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On the Motion of an Ellipsoid in a Fluid. 

By Thomas Craig, Fellow of the Johns Hopkins University. 

Let us suppose a solid body of the form of an ellipsoid whose mass is so 
distributed as to be symmetrical with respect to the three principal planes 
of the body ; that this body be immersed in an infinite mass of frictionless 
incompressible fluid, the whole system being originally at rest, and that the 
solid be acted upon by any set of impulsive forces— or by any impulse — and 
the system then left to itself; to determine the resulting motion of the solid 
body and of the fluid. Concerning the motion of the fluid we can say that it 
will be subject to a velocity potential, since the entire motion of the fluid 
being due to that of the solid — or to the motion of a portion of the bounding 
surface of the fluid — there can be no rotation in any part of the mass. 

We will now proceed to the determination of the velocity potential <p 
satisfying the equation 

V 2 <?> = 0, 

V 2 as usual standing for the operation 

9! , 3 2 , & 

dx 2 "*" dy 2 "■" dz 2 " 

As we limit ourselves to a simply connected space, the function (p will be sin- 
gle valued. 

Let now if, q, £ denote the coordinates of a point with respect to a set of 
rectangular axes fixed in space, and also let x, y, z denote the coordinates of a 
point with respect to a similar set of axes fixed in the body, we have then 

£ — a + a l X "f" a 2J/ + 0*-3 Z ) 

yi = (3 + p 1 x + (3 2 y + (3 3 z, 
f = y + y t x + y%y + y-iZ , 

the twelve quantities a, (3, y being functions of the time and position of the 
body whose geometric meaning is well known. 

If we denote by u, v, w the translation velocities of the origin of x, y, z 
in these respective directions and by p, q, r the components of angular velo- 
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city around these same axes, we have for the velocity of the fluid particles 

relatively to the body 

dx dw 

-dt=£- u - z 2+y r=:v ' 

1. - d -=^-v-xr + z P = v, 

dz dw 

~di = i~ w ~yp + x i = w • 

If we denote by (nx), (ny), (nz) the angles which the outer normal to the 
surface of the body makes with the axes of x, y, z respectively, we have 

o 

~- = (u + zq — yr) cos (nx) + (v -f- xr — zp) cos (ny) + (w + yp — xq) cos (nz) . 

As the fluid is to be at rest at an infinite distance from the body, the first 
derivatives of $ with respect to x, y, z vanish at infinity ; and since y 2 <?> = 

throughout the space filled by the fluid, and the quantities -— , >~ , --^ are 

continuous throughout the same, we have (Kirchhoff, pg. 225) for ty the 
expression 

<p = ufa + vfa + wfa + pfa + qfa + rfa . 
The six functions fa , fa , . . . fa satisfying the equation \7 2 <p z= and also giv- 
ing at the surface of the body 

-~ = cos (nx) , ^ = y cos (uz) — z cos (ny) , 

2. -~ = cos (ny) , ~- 5 =z z cos (ux) — x cos (nz) , 

~ = cos (nz) , ^ 6 = # cos ( u y) —y cos { nx ) > 

From these equations we see that the functions fa, fa, . . . depend simply on 
the form of the surface of the body, and not at all on its motion, $ being a 
linear function of the quantities u, v, . . . r. 

We shall resume now the examination of the function fa observing that, 
in addition to the equations already given to be satisfied at the surface of 
the body and throughout the space filled by the fluid, it must also satisfy 
for X x =z go the equation 

Denoting by Q. the potential of the ellipsoid upon any point x, y, z, we have 
for a point at the surface of the body 

n = const. — 7t (Ax 2 + By 2 + Cz 2 ) 
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(Kirchhoff, pg. 130), where 

4. A = aicf (a ,f m , B = a b e f-K— 1 C=abcf^f^ 

and for the derivatives of II (ib. pg. 225.) 

d 2 ti 

dndx = 2n ( 2 ~ A " > C0S M' 

5. a ~ = 2jt(2-B)cos(^), 

-„-~- = 2n (2 — C) cos (nz) . 

Comparing these expressions with the values previously obtained for -S- &c, 

we have 

_ 1 dti _ 1 dQ _ 1_ 3£ 

b - ^ — 2ff (2 - 4) "a» ' &-2tt(2~ B) dy' ^ ~ 2^{2~^~C) dz " 

All the conditions that $ 6 must satisfy with the exception of 

-^- 6 zz # cos («y) — y cos (war) 

will be satisfied if we write 

when _/V~ 3 is an arbitrary constant ; Kirchhoff shows that this condition will 
also be satisfied by a proper determination of N z ; the value which he obtains 
in a very simple manner is 

a 2 — h 2 

Q TIT U, U n 

3 ~ 2tt [2 (a 2 - 6 2 ) + (A-B) {a 2 + 6 2 )] ' 

If instead of ]Y t , 2f 2 , JY 3 equal to the values here given we write 2rtN u &c, it 
will be a little more convenient, and we now have for q> 
Au Bv , Ow 

+ IT 1 {B— C) yzp + N 2 (C— A) zxa + N Z (A- B) xyr . 
The quantities A, B, Care given by the expressions 

A = — bc-7 X 1 , 
da 

t> d«h 

B ■=. — ca - oT - , 
do 

C — — ao -7T- , 
do 
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where o t = f 4 = , 

J V(a 2 + A)(& 2 + A)(c 2 + A) 

and for convenience of future reference 

C dl 2 

J V (« 2 + ;, 



6)3 



=/: 



V(« 2 + a)(& 2 + a)(c 2 + a.)' 



V(« 2 + 4)(^ + ^)(c 2 + ;,,) 

We can arrive at the final form of <|> in a slightly different manner, which 
will bring in evidence some interesting properties of the quantities that we 
are dealing with. For simplicity, suppose that the ellipsoid moves with a 
constant velocity u in the direction of x, then, as before, 

(?) = u$ t , 
with the prescribed conditions for <p y , will satisfy all the conditions of the 



motion, and will give again 



1 dii 



2k (2 - A) d-v ' 

Write ^.^.4=^. £o> £> 



dk t 



Now A + B + C— abc \ C . 2 , . . . 

\ J (a 2 + A)V(a 2 + A)(6 2 + A)(c 2 + A) 

1 r° d ^ , p Mi \ 

- r J»(b°+ A) V (a 2 + ;:,)(6 2 + A)(c 2 +T) ^ (c 2 + A) V(a 2 + AH* 2 + A)(c 2 + A) J 
_ flJc f <«» r(6 2 + A)(c 2 + A) + (a 2 + A)(c 2 + A) + (« 2 + AX^Wl 

c J V(a 2 + A)(6^ + A)(c 2 + A) L {a 2 + A)(6 2 + A)(c 2 + A) J 

_ 2 . abe 

abc ' 

4 + 11^4 = 1, 

now 

1 J-. 1 3i2 0J0 /•» a;r?A 

* l - 27r"(2 - i) ^ — ~~ (2 - J.) J x V(a 2 + A)''W + Aj(c 2 + A) ' 
where % is the positive root of 

substituting for .4 its value, we have again 

a;c?A 
xy/(a* + A) 3 (6 2 + A)(c 2 + A) _ 1 r<* .w//, 



/: 



f » t/A 2 ~ J, + C, J x y^ + /,y(6 2 + aW+ A) 

J V(« 2 + A) : W + A)( e2 + A) «*« 
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and consequently 

-., , u r°> xdX, 

* ~ B, + Cj x V(a 2 + W(b 2 + I{p + A,) ' 
From this value of ^ we can readily find for u, v, w the values 

u f c" dl l x 2x 



XJ = 



- i r 



12 . A L(a 2 + Z) 2 + (^ + Z) 2 + (« 2 + X) 2 J J' 

_ u x 2z T x 2 y 2 z 2 -j-i 

_ u x % r a; 2 y 2 2? ~j _1 

w - ~ a+c, v"(a 2 + Z f(6 2 TzF+B * oT + * ' iw+xt + (T 2 +^) 2 + w+Vy\ ' 

3y 2a; r x 2 y 2 , z 2 -i-» „ . 

Since ^ = -^ ,-r- rr + rrj-^ rj + FT - ! T 2 ) &C -> &C - 

dx a 2 + x U a + t) ( b + Z) (° + X) J 
If we suppose now that 

$ = l«fc + V<?> 2 + W0j , 

i. e. that the body move in any direction with the component velocities u, v, to, 
we see that all we have to do is to determine separately <2> 2 and <2> 3 and, substi- 
tuting in this expression for $, have the motion completely determined. 
Obviously the functions q> 2 and ^> 3 are given by 

_ i r x ydh 

cp ' 2 ~ a,+ aJ x vv + ;>0(6 2 + xjXc 2 + ; u ) ' 

1 /»» zdki 

^ ~ A, + bJ x \/ (a 2 + ^){b 2 + /^{g 2 + /^f ' 
and thus for a simple motion of translation we have 

u /»« xdly v /»» ye//.! 

14. ~" ^+^' ^ x \ / (« 2 +/ 1 ) : \6 2 + A)(c 2 +/,) Cl + ^ l ^ * V(« 2 + ^)(6 2 +/ 1 ) 3 (c 2 +/ 1 ) 

+ a. + bJx V(f7+ /,)(6 2 + l3p + ^f ' 
Similarly in the case where there is rotation we find the values 

1 1 4 r x d«i 



jV - = 2 , 






Ue5 2 - v c + a V 

at - 1 x r - f" rfgi 



L!a -(* + *>)' 
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So that finally we can write for the general value of $ the following 
u r x xdli v r<* ydX l 

$ — B.+ gJk V(a 2 + ^) 3 (6 2 + ; 4 )(c 2 + XJ + A^+Ci J x V(a 2 + ^)(6 2 + ^) 3 {c 2 + V) 
w /•=» zdXi 
+ 27+ Bj x yV + /,)(6 2 + /,)(c 2 + A) 3 
(6 2 — c 2 ) p /•» yzdXy 

16. ^ +C ^ s ~a6e 



)(6 2 + A)\c 2 + /,) 3 
(c 2 — a 2 ) g /•» zxdki 



{cr — a') q /•» 

^ 4_ ^ jj _ !_ •> x V(« 2 + W + AXc 2 + /0 s 
1 ^ ; 2 aba 



(a — 6 ) r /•=» jm/A^ 



*) 



The approximate values of the quantities .4, 7?, Ccan be obtained by expand- 
ing a>! in an infinite series and then performing the indicated differentiations, 
which is the method employed by Olebsch for determining their values ; we 
can also obtain them directly as elliptic functions capable all of being reduced 
to depend upon the ^-function. If we take %i, % 2 , %% as the amplitudes of 
three elliptic integrals 

17 - f-^L_ 6 - C dy - 2 - f -^ 3 - 

17 - ^~JA(k, x y v *-Jj(k >X2 )' ^-jA{k,^ 

we shall be able to give A, B, C as elliptic functions of any one of these three 
quantities X , 2 , 3 and could determine f» as a function of all three, i. e. 
supposing the proper relations to exist between these quantities and the 
parameters \, X 2 , /l 3 ; these relations must evidently be such as to determine 
0i, $2, #3 as the so-called Lame or ellipsoidal coordinates. First to determine 

A, B and C, write 

h? — x 2 a 2 — c 2 

18. ^i = c 2 — 5— , when h 2 z= — -. — 

' ' XT & 

x having all values between zero and infinity. This gives us 

x 2 + 1 



a 2 + % x ~ (a 2 — c 1 ) 



i 



_ (6 2 -c 2 )x 2 + (a 2 -c 2 ) 
IQ ° ~r A i — 77 j 



'H 








X 2 


-*1 


— 


a 2 ~ 


c 2 




x- 


> 




d\ 


— 


2 


(a 2 - 
x :i 


-c 2 ) 



dx 
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a? — b 2 
Assume -^ — § — ^ 2 an ^, consequently, for the complementary modulus 

b 2 — e 2 . ... 2abc 

:=^ 2 ; making these substitutions and, for brevity, writing i" = 



we have 

J °(x 2 + 1)W 1 + k"a? 
20. s=I j7 Xl ' {X 



o (1 + k' 2 x 2 ) V^ 2 + 1)(1 + /i'V) 

c=// o " "**" 



V(a* + l)(l + A'V)' 



7Z 

Again, writing ^^tan;^, where Xi ^' ies between and -^ , we have 



readily 

»• s =if: U«b)- 1 )jws= / /K*k- 1 ) *. 

r | sin 2 y t (/ft _ r | sn 2 ^ 
Jo cos 2 ft J (A, ft) Jo cn 2 ^ Wl ' 
Now we know (Cayley, Elliptic Func, p. 15) that 

/Ppm 2 d 1 dd 1 = (l - ^) Q l - Z6i , 
and so we have readily for A x 

22. A= 4{(l-f)«i-^} 

_ J ) *fl t 8.E y(fl t )i 

and also by not difficult reductions 

_ i_ ] e ± d % k \_ &' (e, + k) 

"*' 7,2 1 7. CIA. ~T* J./2 



A 3* ' k' 2 6(di + K) 

o=^{|; log JA+i)-|fc}. 

The expression for C could be given as depending upon the S-function by the 

formula 

H (6, + K)= «r= (ir '- 2 ^ (0 X + JT ) , 

but nothing would be gained thereby. We may just notice the forms which 

immediately present themselves for the constants A % , B 2 and C 2 , but will not 
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2 
integrate the equations. For brevity, we will here denote - — , by /', and 

(af — c^ 
we have by the same transformations employed in reducing the quantities 
A , B and C , 

a — j< c x ! dx r, r si " 4 My 

2 ~ J (1 + ^V(S +T)(1 + FV) ~ J cos 2 ^ zJ 3 (A, Zl ) ' 

23 B =1'C xV/a; = I'C — 7y l/l 

2 J l* 2 + 1) V(* 3 + 1)(1 + h n x z ) J cos 2 7! J (k, Xl ) ' 

/> _ jC _aMc _ p rsn^x l dy l 

° 2 ~ J(^+l)(l + A'Vv( !B *'+lXl'+*' s i 2 )'~ •> J ^Z0" ; 
and these, by virtue of the equation ^ rz amf^ become 

24 ^= / /d ! M;'».' ^ =/ .r s " ! ">3,*- ^j^d:*- 

The values of x, y , and z can, of course, be expressed by means of the func- 
tions dx, 2 , 9-z . The values of x 1 , y 1 and 2 2 , expressed in terms of the para- 
meters %xi ^2) ^-3! are (see equations 50. in article on "the Motion of a Solid in a 
Fluid; " this Journal, Vol. II, JN T o. 2.) 



_ (« 2 + ;, )(a 2 + A 2 )(« 2 -M,,) 



(a 2 -6 2 )(a 2 -c 2 

o* ._ (&l+ji)(y+j »)(y+j.) 

^°' y - (6 2 -c 2 )(& 2 -a 2 ) 

(c 2 + ;,)(c 2 + / 2 )(c 2 + i 3 ). 



Z 2 = 



(a 2 -o 2 )(6 2 -c 2 ) ' 

performing here operations similar to those that have already been employed 
in transforming to dx, we have, since 6 2 and 6 S sustain the same kind of 
a relation to /l 2 and a, 3 that Oj does to /t t , 

* 2 = (a 2 - , 2 ) ^^ , X = V«^ M.™& ; 

26. f = (a*~c>) d ^ C ^^ 26s , y=Va^? dn6 ^^ ; 

, cn%s4dn ! & /.j , cnflj sn0 2 dn0 3 

Z 2 = (a' — c-) - 2 7i > 2= V«-— 6'" - „ — „ 

With these values of x,y, z, A, B, C, A.,, B. n G 2 , we could now 
transform ty so that it should be given as a function of elliptic functions, but 
nothing of interest could come from that operation in the general case where 
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none of the velocities vanish, so we will not attempt it. For the case of Jc zz , 
or an oblate ellipsoid we have readily 

27> A = B = ^-($ 1 - sin 6, cos 0,) 

C= I (tan 6,-6,). 
Suppose Jc z= 1 , i. e. a prolate ellipsoid, then we have 

%i = gdOi 
and A (Jc, %,) = cos %, = cg0j , &c. We have in this case, as is well known, 

X = log tan (— ?t + -y £i) ~ log tan (y »* + y gd^i) , 
and, passing to exponentials and reducing, 

-* - °g^» 

and thence 

e" 1 — f-* 1 

S g^l = 0. I 9. = tanh ^11 

c § -0 i = 6 '. + 6 -'i sech 6l • 
Substituting in the general values for A , B and C and we have 

28. A - lfsg'6, d6, ~ if (I — cg 2 ^) tf0 x , 

or, 

4 = 1(6,- sg-00 =z 7, (0 X - tanh X ) ; 
and similarly, 

= /,/sinM), «» = -i { ^ + (1-29.) j. , 

I, being what 7 becomes for b zz c . 

The motion of the fluid particles relatively to the body is given in recti- 
linear rectangular coordinates by equations 1. , but, as for the problem under 
consideration it is more convenient to use elliptic coordinates, we will con- 
sider equations 39. of the article on " the Motion of a Solid in a Fluid" above 
cited. 

We will first assume the case where the ellipsoid has merely a motion of 
translation in the direction of one of its axes — the axis of x for example ; for 
this case we have 

vzzwzzp = 2 = t ■=. 0, 
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and the resulting equations of motion, 

™ dh _ df dx 

* dt~dk u 3V 
&2 dx, _ d^ _ u dx _ 

dt dh dh ' 

These are the equations- given and discussed by Clebsch. We have also for <p 

A __ Ax da 

$-A=2* U -A^2'dt> 
and 

dx _ 1 x dx 1 x dx 1 x 

ai; ~ t o» + *i ' ax 2 ~~ t o ! "+t 2 ' ai 3 ~~ t a 2 + 4 ' 

these can cleai'ly be expressed simply in terms of the elliptic functions^ — but 
nothing would be gained by it, as the resulting equations of motion will give 
rise to integrals of a higher order than elliptic. Making these substitu- 
tions in the differential equations of motion and they become after simple 
reductions 

(h-kWi-h) d^_jx_ 2 f 2 \ 

(a 2 +A 1 )(6 2 +A 1 )(c 2 +A 1 ) ' da ~ a' +X, ' A-2 \ {A-2) , v /(a 2 +/l 1 )(6 2 +/i,)(c 2 +/i 1 )J ' 
(4 -^3)^2-^1) dh_ 2x 



30 



Making 



(a 2 +X 2 )(b 2 +X 2 ){c*+X 2 ) * da ~ a 2 +X 2 ' A-2 ' 

(h ~h) {h -k) dX 3 _ 2x 2 

(a 2 +A 3 )(6 2 +A 3 )(e 2 +A 3 ) * 3« _ a? +k ' A-2 ' 

Ax da 



A-2 (X.-X^-X^X.-X,) 
we can write these equations in the simpler form 

(a, - a 3 ) d9 = (62 + Xi) p + h) • - 



= <?$ 



{A - 2)V(a 2 + A,)(6 2 + AOlo 2 + A,) 
31. (X, - a,,) dV = (62 + j^p + 4) , 

(a* - ^) rf* = ( jr+j^^qr^ ' 
Multiply these equations respectively by 

v -|- Ai , C "T" A2 ) C 2 -p A 3 , 

and again by 

J 2 + A ls & 2 + A 2 , S 2 + A 3 , 
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and add each of the sets thus obtained, observing that the quantity 

1 



can be written 



(A - 2)V(a 2 + XJ(b* + \){<? + A,) 



1- 



1 



(4-2)V(a*+W + W + *i) 



and we have for the equations of the motion of the fluid particle relatively to 
the body 

y — 



b 2 +?n 



32. 



since 



9 dz dX t 

T _ 7 2 TXi 



(4-2)V(a» + ^(6* + A I )( * + ^) 



1 - y (^ - 2 V(a 2 + ^)(6 2 "+ ^)(c 2 + ;,) 



* 2/ ~6 2 + A 1 ~^ 2 + ^ 



+ « 



<H* 



b 2 + I 



, &C. 



From these we see that the path of the particle lies always in a plane passing 
through the axis of x ; by taking values of \ so large that the fourth powers 

la? IP Ic 2 

of the ratios \'— , \r- , \ T , can be neglected, Clebsch shewed that the path 

h A A x 

of particles very rapidly approached straight lines. By subtracting the 
second of these equations from the first and integrating we have 

33. l g^-=-(5 2 - C 2 )f - ^ r. 

J (6 2 + ioco 2 + k) [ i - Y {a - 2)v[^w+w~-no } 

This can be expressed in terms of the elliptic functions, for we have 

r- dh _ 1 dQ l 

J V(a 2 



V(a 2 + W + ^)(c 2 +; 4 ) 

1 dOi snU 



2 Va 2 



d\ = 



2 Va 2 -c 2 'cnftdnfo' 



(& 2 + ^)(c 2 + ^)=(« 2 - C 2 ) 2 . 



cnWft 
sn 4 0, 



Making these substitutions, we have 

sn 7 t dft 



34 log(-^) (a2 - c2)i = T 



cnfydn'Gi 



1 



. k 2 {a 2 -c 2 )i 



■Zdi 



-1 



/„2 ^fcn&dnft 

\ a ~ c ) Tr\ 
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a non-integrable expression. Nothing of interest could be obtained by a fur- 
ther examination of this case, and so we shall leave it, calling attention, how- 
ever, to the fact that the coordinates of any fluid particle are expressed as 
functions of the parameter of the confocal ellipsoid upon which it lies. The 
same is also true in the case of a simple rotation around one axis — as the axis 
of x\ the resulting integral equation of the path of the particle comes in this 
case also in a form that is too complicated to admit of any discussion. The 
motion of the particles, very near the surface of the body, would be, in an 
ordinary liquid, a very important point to determine- — but, as in that case, the 
equations will present still greater difficulties than they do in the case of a 
perfect fluid, it does not seem as if that problem is susceptible of solution. 

We shall now take up the second and more interesting part of the gene- 
ral problem and investigate the motion of the solid body. As there are no 
external forces acting upon either the body or the fluid, the motion of the 
system is entirely due to that of the body, and, in consequence, the total 
energy of the motion is a homogeneous quadratic function of the quantities 
it, v, w,p, q, r with constant coefficients. On account of the assumed sym- 
metry of figure and distribution of mass of the body, i. e. with respect to 
three mutually perpendicular planes, the expression for the energy of the 
system will contain only the squares of it, v, &c, and may be written 

2 T= a n ii 2 -f- a 22 v 2 + %jW 2 + #14? + a*, t + ««c^ 2 > 
where the constant coefficients a, 5 depend upon the density of the fluid, the 
form of the surface of the body, and the moments of inertia about the axes. 

Write for brevity 

dT dT 

— — U — - — R 

ou or 

The Kirchhoffian equations for the motion of the body are now 
~ = q W- r V, 3 d f = v W- iv V+qIi~ rQ , 

35. ~) --rU-pW, d -^ = wU-itW+rP-jpB, 

^- u V-vU-\-pQ-qP. 

From these equations we have the well known integrals 

2T= const. , 

36. U 2 + V 2 + W 2 - const. = L , 

UP + VQ +WP = const. = 21 . 



dv 
dt ~ 


zrU- 


-pw, 


dW 

dt ~ 


-fV- 


-qJJ, 
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For the assumed case the equations of motion become 

«n -Z = a mm — <hW , a u J^— (a 33 — a 22 ) vw + {a m — a 55 ) qr , 

37. a 22 ~ = a u tw — « 33 wp , « 55 ^ = (« u — <? 33 ) wu + (« 44 — a 66 ) pr , 

« 33 -^ = a-22»p — a n uq , a<x^ t = O22 — a n ) uv + (a 55 — a u ) pq_ . 

Before going on to the integration of these equations for special cases, it 
will be desirable to determine the values of the constants a i5 . These quanti- 
ties are made up of two parts depending, as they do, upon the energies of the 
solid and the fluid ; that part of a n which depends upon the solid merely is 
m, or the mass of the solid; and the same is true for a 22 and « 83 : the parts 
of « 44 , a m , a m depending upon the body are 

f(& 2 +c 2 ), f (c*+a>), f («■ + **). 

For the parts depending upon the energy of the fluid we have, denoting any 
one of them by a' tj , 

a' tj = pffy ^ da 

(see article above cited, or Kirchhoff 's Math. Pkys.) ; the integration is of 
course to be extended over the entire surface of the ellipsoid. 

The values of these constants have been already obtained by several 
writers, and I will merely write them down — they are : 
, _ p A__ 

38. a'22 = m — 77-j — r , Po denoting the density of the solid. 

Po + A 

P o 



p A + B ' 



'° 2-C+f-^{A + 2+Q 



39. ^f^+^jl + A. °- A 



C-A+^(A + B+G)[' 



, m , 2 7 ox (.. . p A — B 



A-B + -^ 2 {A + B+0 
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The first three of these are the same in form as the values deduced by Mr. Fer- 
rers in a recent number of the Quart. Jour. ; the last three may be obtained 
from his values or from the values already given for N x , N 2 , JV" 3 . 

The values of these constants are thus seen to depend upon elliptic func- 
tions, and they can be readily determined from the previously given values 
of A 2 , B 2 , C % , or they may be reduced to depend merely upon A, B, C. Con- 
sider for a moment the last three, viz : a u , a 55 , a m ; they are the moments of 
inertia of the ellipsoid, with respect to its axes, increased by certain quanti- 
ties depending upon the density of the fluid and also upon the density and 
ellipsoidal shape of the immersed body. Conceive another ellipsoid whose 
semi-axes are a x , b x , c x and of mass m; the magnitude of the axes to be 
determined from the conditions 

•y (6? + cf) = «* , 

40. ^{<* + af) = a B , 

~ (a\ + b\) = a m . 

Denote for brevity the fractional quantities in the expressions for a u , a 55 , a K 
by £i > £2 > £3 > then we have for the determination of a\, b\, c\ the equations 

bl + c\=(l/+c')(l + Q, 

41. c ? + a? = ( C 3 + a 8 )(l + f 2 ), 

a\ + b\ = (a* + F){l + Q; 

these give readily 

„ 2 _ „* , (« 2 + 5 2 )c3 + ( c 2 + o') c.-(&' + ^:i _ . , „. 

«! _ a -\ -- =: a -j-ti, 

42. Sf-S 2 j (« 2 + ^ 2 )C3-(o 2 + a 2 )C 2 + (^ + o 2 )C 1 ---y | 7 2 > 

^_, 2 , -(» 2 +^ 2 )C3+(c 2 +a 2 )c 2 +(5 2 +c 2 ): 1 _ . 

and we have also for the new density 

( abc _ abc 

If the ellipsoid is oblate 

a\ + b\ = a 2 + b- , 

if it is prolate, 

b\ + c\ = ¥-\-c\ 
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A further study of the properties of this ellipsoid would undoubtedly prove 
of interest, as it seems to possess, as will be seen shortly, properties some- 
what analogous to " Poinsot's ellipsoid " in rigid dynamics. 

The transformation to Lagrangian, or generalised coordinates, gives rise 
to several interesting results, and before going further with the discussion in 
Eulerian coordinates it will be desirable to make this transformation. I will 
write down at once the values of u, v, w, p, q, r in terms of the generalised 
coordinates x, y, z, 6, <£>, 4"> 6, $, 4>, as they are obtained in Thomson & Tait's 
Nat. Phil., (new edition), and also in " The Application of Generalised coordi- 
nates to the Kinetics of a Material System" by Watson and Burbury ; they 
are 

u = (cos d cos q> cos 4' — sin 6 sin 40 x 

+ (cos cos <p sin ^ -|- sin <p cos 4) y — sin 6 cos <pz , 
a. v = — (cos 6 sin ^> cos 4- + cos <p sin 4^) x 

+ (cos 6 sin <£ sin 4- — cos ^ cos 4-) y + sin 6 sin fz , 
w = sin 6 cos ^x + sin sin 4# + cos 6z , 
and for the rotations the well known values 

p =z sin q>Q — sin d cos $4- , 
,3. q — cos cph + sin sin ^ , 

r = <?> + cos 6^ • 
The Lagrangian equations of motion are also, no external forces acting, 

ara*- u ' a< a#- u ' ar~ai- u > 

dtdo dd~ ' ar a? a f ' ar ^, ~ ^ - ° ; 

these equations can be found in either of the above mentioned books, or in 
Routh's Rigid Dynamics, and are too well known to need any explanation. 
The values of x, y, z could be obtained by direct solutions of equations a, but 
that would be a long process, and it is obvious from their forms that the 
desired quantities can be otherwise obtained. Multiply u by cos 6 cos <p, v by 
cos sin <?>, w by sin 6 ; and again multiply u by sin f, v by cos $ ; and finally 
u by — sin cos <p, v by sin 6 cos $, and w by cos 6 and in each case add the 
results ; we have then 

x cos 4- + y sin 4' — u cos $ cos <£> — v cos 6 sin $ -f w sin 6 
— x sin 4- + y cos 4- = u sin $ -f « cos 4) , 

2 = — u sin cos <?> -f « sin cos <p + w cos 

= (« — «) sin cos <p -f- to cos , 
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the values of x, y and z obtained from these equations are readily found to he 
x = (cos cos ^> cos 4- — sin <p sin 4) u 

— (cos sin $ cos 4 + cos <p sin 4) v + sin cos 4-w , 

44. y = (cos cos <|) sin 4 + sin q> cos 4) « 

— (cos sin <p sin 4 — cos <p cos 4) f + sin sin 4-w , 
z — — w sin cos (|> + ^ sin cos $ + w cos . 

Having thus determined x, y, z as functions of 0, <p, 4, if we can now, from 
equations /?, determine these latter quantities as functions of 0, q>, 4>, we shall be 
in condition to integrate these differential equations and so determine the path 
described by any point of the body during its motion through the fluid. We 
know that whatever motion the body may have at any time that we can 
conceive the motion as clue to a properly applied impulse at the beginning of 
the motion. It will simplify the work to assume the axis of the impulse as 
coinciding with one of the axes of reference, as the axis of z ; then calling L, 
see equations 36, the total momentum of the impulse we will have from equa- 
tions a — the components of momentum in the directions of x and y being zero, 

45. — — , — -r =z , -t = L, (a constant) . 

dx dy dz 

We have further, Thomson & Tait, § 221, for the Eulerian components of 

momentum a n u, a 22 v, a 3i iv 

a n u cos $ — a 22 v sin ty — — L sin , 

46. a u u sin $ + a 22 v cos <p = , 

a 33 w — L cos ; 

from these follow 

a n u ■=. — L sin cos <|> , 

47. a 22 v = L sin sin $ , 

a S3 w = L cos . 

From these we have, for the first three terms in the expression for the energy, 

a n u 2 + «22^ 2 + <%w 2 = 

48. D Tsin 2 ( C0 ^ + ****-) + ^1 = G . 

L \ a n a- n J o 33 J 

For the last three terms containing the squares of the angular velocities, 

49. a u (0 sin $—4 sin cos <|>) 2 + a- o5 (0 cos $ + 4 sin sin <p)' 2 + « 66 (4 cos + $} 2 ; 

this quantity is identical with the expression for the energy of a free rigid 
body rotating about its center of inertia (Thomson & Tait, p. 314), the 
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moments of inertia being a u , a 5S , « 66 . We have for the total kinetic energy, 
m 2T= G -f a u (6 sin <p — 4- sin 6 cos $) 2 -f- « 55 (0 cos <j> + 4 sin sin $) 2 

+ « 66 (4cos0 + 4>) 2 ; 
from this it is clear that the motion of the ellipsoid in the fluid is identical 
with the motion of another ellipsoid of semi-axes a lt b 1} c x rotating about its 
center under a potential 6r, or, the entire kinetic energy of the system of 
fluid and ellipsoid is equal to the entire energy of an ellipsoid of semi-axes 
«! , b x , <?i and density p' rotating about its center under a potential Gr. 

If the ellipsoid be of revolution around the axis of z the expressions for 
x , y , z become much simplified ; we have in that case 

"ll— #22 ) #44 — #55 ) 

and, as we have already seen, 

«6G = -y (# 2 + & 2 ) • 
The values of x , y , z are now 

x zz — i sin 6 cos cos 4 ( ) , 

51. y z= L sin = cos 6 sin 4- ( ), 

Vein a i3 / 

z- L sin 2 e(— -—)+ — . 

[A very interesting problem is solved by Thomson & Tait, viz : to find 
the motion of a solid of revolution in a fluid which so moves that the axis of 
revolution shall always be in one plane — the body has a motion of transla- 
tion in the direction of the axis of x and an angular motion around this axis — 
and it is found that the resulting equation of motion is the same in form as 
that of the common pendulum, the mass, moment of inertia and length of the 
pendulum depending upon the mass and figure of the body and the density 
of the liquid.] 

The substitution of the relations a n z= a 22 and a u z= a 55 in the expression 
for 2T cause this quantity to become 

52. 2T- L 2 sin 2 6 (———)+—+ a u (0 2 + sin 2 eft) + a m (<?> + cos 4) 2 

containing neither <p nor ^ , and, consequently, cause the fifth and sixth of 
equations y to become 

or, dT . dT 

' — = const. — =z const. ; 
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equating then to constants these derivatives of Tand we have 
53 <?> + cos Oi = fr , 

a u sin 2 64> + a e ^i 1 cos = [i 2 » 
,«! denoting the angular velocity around the axis of z and ft 2 the component 
angular momentum about the same axis. The equation 

d_ dT_dT 

3< ' 30 ~ 90 ' 
now merely expi'esses what we already know, that T, the kinetic energy, is 
constant. 

The expression for the energy may now be written 

j} sin2 6 (-1 - - 1 -) + ^ + a u (0 2 + sin 2 Qp) + a, 4 t\ . 

V*ll "S3 / a 33 

The steady motion of the solid is easily obtained from the original expres- 
sion for the energy, together with the values already given for u, v, w, p, q, r 

■ ■ ■ " ' dT 

and for x , y and z . We must have 0=0 and consequently ^-- = ; since 

dT dT . 

7T- and — are identically zero, they of course do not come into the question, 

we must also have 4 1 = const., say v; . 
d_T_d_T du dTdr_, 

dQ~du'dd + drdd' 

dT__ dT_ 

du - °' nU ' dp - a ^ ' since a n = a S3 

dT dT anda 55 = a 60 . 

dv dq 

dT dT 

di 
du 

m 

du 

di 

dw 

36 
dp 



dv - ° wV ' dq ~ a5 '° q ' 
dT 

dw or 

■=. — (sin cos <£> cos ^x + sin cos <p sin ty + cos cos <pz) , 

du 

= (sin sin $ cos 4* + sin sin 4) sin tyy + cos sin q>z), 

06 

— (cos cos ^ + cos sin 4^ = sin 0z) , 



z= — cos cos 4>^ , 
9p 

5 ? = cos sin ^ , 
dp 

9/- . „ : 

d6=-* mH > 



71 
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dT 
Substituting all these values in ^ , and also giving p and q their values in 

06 

the generalized coordinates, and this becomes 

dT . . • 

■^r = — a n (u cos <p — v sin <j!>) (sin cos 4# + sin ty + cos 6z) 

of) 

+ a S3 w (u cos <£> — v sin <?>) 

— a^ (sin $0 — sin cos <?>4) cos cos <|>4 

+ a 44 (cos <pd m sin sin <p4>) cos sin $4 

+ a m r sin 04 , 

since 

u cos <?> — v sin <£> = # cos 4 cos + y sin 4- cos — z sin 

= ( ff 33 — ^n) ( M cos <?> — w sin <£>) w 

+ a 44 sin 2 04 2 — <W cos 04 

= D sin cos ( ) + « 44 sin cos 6^ — a^r sin 04 . 

c \ct33 ctu / 

Calling now the integral of 6' = , 6 = 6', then cos rr cos 0', we had 4 = v > 
therefore, cos 4 = cos qt, &c. The equation ^- = , becomes now 

54. « 44 cos 0V — a m ryi -f Z 2 ( ) cos 0' = . 

which gives 

1} (L _ L\ cos 0' 

\«33 _ a ll/ 

the value of the angular velocity necessary to the maintaining of steady 
motion when the body is projected in the direction of its axis. The equa- 
tions for x, y, z now give after integration 

x = — ( ) sin 0' cos 0' sin nt , 

q \0r33 «„/ 

56. yr — ( ) sin 0' cos 0' cos >?£, 

^ \<*33 °11 / 

\ <% a,, / 

or 



55. r = ^» a '^ + ^ cos 0>7 , 



e;7 T.~ — ^- an *?^ 



a; 
2/ 



a = L cos 2 0' f — - --) t + — , 

W Oil/ Oil 
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the equations of a helix whose projection on the plane of xy is the circle 

58. x 2 + f = K (- -) 'sin* 0' cos 2 9. 

rf Vaits a n J 

This is in accordance with what was proved in the paper already referred 
to, viz : that a body moving in a fluid has an infinite number of possible 
steady motions, each of which consists of a twist about a certain screw. 

The problem of the motion of a solid of revolution in a fluid has been 
very fully solved by Kirchhoff in Borchardt's Journal, Bd. 71, and by others 
since that time. I append a list of the more important articles on this sub- 
ject : 

Olebsch : Die Bewegung eines Ellipsoids in einer Fliissigkeit. Orelle's Jour- 
nal, Bd. 52 and 53. 

Kirchhoff: Ueber die Bewegung eines Rotations Korpers in einer Fliissigkeit. 
Borchardt's Journal, Bd. 71. 

Ferrers : The motion of an infinite mass of water about a moving Ellipsoid. 
Quart. Jour., No. 52. 1875. 

Kopcke: Zur Discussion der Bewegung eines Rotations Korpers in einer 
Fliissigkeit. Math. Annalen., Vol. 12, p. 387. 

Weber: Anwendung der Thetafunctionen zweier Veriinderlichen auf die 
Theorie der Bewegung eines festen Korpers in einer Fliissigkeit. 
Math. Annalen., Vol. 14, p. 173. 
The motion of a sphere in a fluid was treated by Dirichlet, and was 

the first article of note which treated of the motion in a fluid of a body of 

given form. 

Washington, July .5, WO. 



